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LOCATION METHOD FOR MOBILE NETWORKS 



Field of the invention 




The present invention relates generally to location techniques. 
5 More specifically, the present invention relates to determination of the geo- 
graphical location of a mobile (i.e. a mobile terminal) within a mobile network. 



10 eled the development of location determining techniques in mobile networks. 
First, different authorities set requirements for the location determination of 
mobile terminals. It is highly desirable that certain authorities, such as emer- 
gency call centers, can locate the calling party as accurately as possible. In 
many countries legislation sets requirements for such location methods. For 

15 example, in the USA mobile location to an accuracy of 50 meters for 67 per- 
cent of calls and 150 meters for 95 percent of calls will be mandatory in the 
near future for handset-based solutions. Second, many of the future services 
provided in mobile networks will be such that they require information about 
the current geographical location of the mobile terminal. 

20 The greater is the accuracy of a location method, the better it can 

serve the application utilizing the location information. This applies especially 
to densely built urban areas where long-range visibility is not possible. The 
accuracy of the location methods is dependent on many different factors, 
such as radio propagation effects. This is a major reason why dense urban 

25 areas form a difficult and challenging environment in terms of location accu- 
racy, with severe multipath characteristics and signal propagating with less 
attenuation along street canyons than through buildings. 



obtained from the network. The measurement results are typically not re- 
30 ceived as coordinates but as parameters which define a certain area within 
which the mobile is located with a certain probability. As the network normally 
provides several parameters for a single measurement, the parameters relat- 
ing to a single measurement are in this context denoted as a parameter set. 
On the basis of the parameter set, geometrical limitations can be determined 
35 which define the geographical boundaries of one or more areas where the 
mobile is most likely located. As discussed below, the determination of the 
location in current cellular networks is widely based on a parameter set in- 



Background of the invention 

There are two major reasons that have given motivation and fu- 



When a mobile is located, successive measurement results are 
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eluding the cell ID and the Timing Advance value, which indicate how far 
from a certain base station the mobile most probably is located. 

The current location methods employ the parameter sets obtained 
by first defining the said one or more areas and then calculating an estimate 
5 for the location of the mobile. History data, i.e. previous estimates, can then 
be used to make the current estimate more accurate. One known method, 
which is based on a least squares scheme, determines a path for which the 
sum of the squares of the deviations from the estimated location points ob- 
tains a minimum value. 

10 A major drawback of the current methods is that they do not pre- 

serve the original data revealed by the parameters obtained from the mobile 
network as to the geographical distribution of the location of the mobile. In- 
stead, with current methods the estimation process moves in an early phase 
to a point-based approach where the information concerning the geographi- 

15 cal distribution is reduced to estimated location points and the final estimate 
is calculated on the basis of these points. Once the process moves from the 
geometrical data to point-based information, it can no longer regain the in- 
formation about the geographical distribution of the location which was in- 
cluded in the parameter sets received from the mobile network. 

20 The objective of the invention is to eliminate the drawback de- 

scribed above and to bring about a solution which enables the accurate loca- 
tion of the mobile through efficient usage of the information obtained from the 
mobile network. 

25 Summary of the Invention 

The objective of the present invention is to find a solution to im- 
prove the accuracy of the current location methods. Furthermore, the objec- 
tive is to achieve a solution not requiring extensive prior labor for achieving 
the improved accuracy, such as field measurements or data collection. 

30 These objectives are achieved with the solution defined in the inde- 

pendent patent claims. 

The invention utilizes location-dependent parameters available from 
a mobile network for determining the location of the mobile terminal. As men- 
tioned above, the parameters relating to a single measurement are in this 

35 context denoted as a parameter set. These parameters indicate the geo- 
graphical areas within which the mobile most probably is located, i.e. how the 
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location of the mobile is geographically distributed over the area of the sys- 
tem. 

In order to improve the accuracy of the system, the invention util- 
izes the geographical information provided by the parameter sets by forming 
5 matrices which indicate the probability distribution of the location of the mo- 
bile. A matrix is formed for a parameter set received from the mobile network. 
Each element of the matrix corresponds to a certain geographical area and 
contains a value which indicates the probability that the respective mobile is 
within the said area. At least one matrix formed for a mobile is stored as 

10 history data to be used in connection with a subsequent parameter set re- 
ceived for the mobile. When such a parameter set is received, the matrix 
corresponding to the current parameter set is formed or retrieved, and the 
elements of the matrix stored earlier are updated. In the updating process the 
effect of the movement of the mobile on the matrix elements is taken into 

15 account, the movement being the estimated movement occurring between 
the measuring moments corresponding to the respective parameter sets. The 
location estimate is then determined on the basis of the element values of the 
matrix corresponding to the current parameter set and the element values of 
the matrix with the updated values. 

20 Thus the invention utilizes history data in the form of at least one 

previously formed matrix, which is updated according to the estimated 
movement of the mobile. By employing the matrices, the original data re- 
ceived from the mobile network on the geographical distribution of the loca- 
tion of the mobile can be retained. The accuracy of the system can hereby be 

25 improved by processing the current matrix by one or more previous matrices 
associated with the mobile in question. 

In one preferred embodiment of the invention, map information is 
employed in determining the element values of the current matrix. The map 
information preferably includes information about the type of surface in the 

30 geographical area related to an element. The description can be simple, such 
as road, forest, river, etc. The element values of the current matrix are then 
weighted according to the type of surface in question. 

In another preferred embodiment map information is employed in 
updating the element values of the matrix formed in connection with a previ- 

35 ous parameter set. 

In addition to the improved accuracy of the system, a further advan- 
tage of the invention is that it does not require any special mechanism or 
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laborious steps for forming the matrices. The matrices can be formed "on- 
the-fly" based on the parameter set received, or they can be formed in ad-, 
vance and stored in a database. In the latter case the parameter set received 
is used as a search key for retrieving the correct matrix from the database. 

A still further advantage of the invention is that it is not dependent 
on the network implementation but can be applied to any network where at 
least one parameter dependent on the location of the mobile is available 
enabling a matrix to be formed. The method can therefore be used on top of 
network implementations based on different location-dependent parameters. 



Brief Description of the Drawings 

In the following, the invention and its preferred embodiments are 
described more closely with reference to the examples shown in Figures 1 to 
9 in the accompanying drawings, wherein: 
15 FIG. 1 illustrates the location-dependent information available from a 
typical cellular network utilizing omni cells, 
FIG. 2 illustrates the location-dependent information available from a 

typical cellular network utilizing sectored cells, 
FIG. 3 illustrates the system of the present invention, 
20 FIG. 4 illustrates a sample matrix formed on the basis of a parameter set 
received from the mobile network, 
FIG. 5 is a flow diagram illustrating the method of the invention, 
FIG. 6 is a flow diagram illustrating the processing of a preceding matrix, 
FIG. 7 illustrates how the matrix of FIG. 4 is changed when it is updated 
25 once, assuming that all directions of mobile movement are equally 

probable, 

FIG. 8 illustrates how the matrix of FIG. 4 is changed when it is updated 
once, assuming that the mobile is traveling directly eastwards, and 
FIG. 9 illustrates one way of forming the matrix which corresponds to a 
30 parameter set received. 

Detailed Description of the Invention 

As mentioned above, the method of the invention applies to vari- 
ous kinds of location-dependent information. Depending on the particular 
35 cellular system, the location-dependent information provided by the network 
can be signal strength or signal delay, for example. The determination of the 
location in current cellular networks is widely based on the Timing Advance 
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value, because the Timing Advance value is directly available from the net- 
work. Therefore, Timing Advance is in this context used as an example of the 
location-dependent signal information available from the mobile network for 
location determination. 
5 As is known, Timing Advance indicates how far the mobile most 

probably is from the base station. FIG. 1 illustrates the location dependent 
information provided by a network with omni-directional base station anten- 
nas, whereas FIG. 2 illustrates the same in connection with sectored cell 
sites. The network typically provides the Timing Advance information as the 

10 minimum and maximum distance from the antenna (R min and R max ), in which 
case the mobile terminal is with a certain probability between these limits, i.e. 
the hatched area A in the figures forms the Timing Advance zone defined by 
said limits. In addition to the Timing Advance information, the network pro- 
vides the cell identifier CID, which identifies the cell in which the mobile ter- 

15 minal is located. This information can be given as the coordinates of the cell 
site. The network further provides an identifier for identifying the mobile in 
question from among the other mobiles and a time stamp indicating the mo- 
ment of location measurement. In the case of a sectored cell, the network 
also provides the sector information. Thus for each location determination the 

20 network provides a parameter set which commonly includes the following 
information: cell identifying data, such as the coordinates of the Base Trans- 
ceiver Station, Timing Advance information, such as R max and R m i n , an identi- 
fier for identifying the mobile in question from among the other mobiles, a 
time stamp indicating the moment of the location measurement, and possibly 

25 the sector information. 

FIG. 3 illustrates the key elements of the system according to the 
present invention. It is assumed here that the mobile network is a GSM Public 
Land Mobile Network. Communication between the network and a mobile ter- 
minal MS in a cell takes place via a radio path by way of a Base Transceiver 

30 Station (BTS) 31 . The Base Transceiver Stations are connected to Base Sta- 
tion Controllers (BSC) 32. Several Base Transceiver Stations are usually under 
the control of one BSC, and several Base Station Controllers are connected to 
one Mobile Switching Centre (MSC) 33, which carries out the main switching 
functions of the mobile network. In addition, the MSC connects the mobile 

35 network with external networks. For positioning purposes, the MSC is con- 
nected to a Gateway Mobile Location Center (GMLC) 34, which collects mo- 
bile positioning information in a positioning database 35. As mentioned 
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above, it is assumed in this context that the GMLC supports Cell ID and 
Timing Advance. This means that for each location determination the GMLC. 
provides a parameter set, including the. parameters discussed above in con- 
nection with FIG. 1 and 2. 
5 The parameter sets available from the mobile network are proc- 

essed in an accuracy server 38, which receives location requests from exter- 
nal objects, such as service applications residing in the network to which the 
accuracy server is connected. In order to be able to determine a location 
estimate on the basis of a parameter set received, the accuracy server forms 

10 a matrix on the basis of the parameter set or retrieves a matrix corresponding 
to the parameter set from among a plurality of matrices formed in advance 
and stored in a matrix database 39a. The accuracy server further uses mo- 
bile-specific history data stored in a history database 39b. This data com- 
prises matrices obtained in connection with the parameter sets received 

15 earlier for the mobile in question. The content of these matrices is discussed 
below. In determining the element values of the matrices, the accuracy 
server preferably further employs map information stored in a map database 
39c, although the use of map information is not necessary. The same applies 
to the use of in-advance calculated matrices, i.e. the matrix database is not 

20 necessary, but the matrix corresponding to the current parameter set can be 
calculated in real-time. 

The information required by the accuracy server for processing the 
parameter sets can also be stored in a single database. The information is 
preferably stored in connection with the accuracy server, although it can be 

25 distributed in the network to which the accuracy server is connected. 

Upon receiving a parameter set, the accuracy server calculates a 
matrix corresponding to this parameter set. On the other hand, if the matrices 
have been calculated in advance, the accuracy server retrieves said matrix 
from the database 39a, and the parameter set is used as a search key for 

30 finding the correct matrix. The matrix contains nxn elements, each covering a 
certain geographical area and each including a value indicating a probability 
that the mobile is within the area covered by the element. FIG. 4 shows an 
example of a matrix comprising 10x10 elements. In practice, the number of 
elements is much greater, e.g. 40x40 elements. Each element corresponds 

35 to a coordinate pair Xi/Yi (i=0...9), which defines the geographical area cov- 
ered by the element. In practice, the size of each element can be 25x25 m, 
for example, in which case a 40x40 matrix would cover an area of one 
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square kilometer. In the example of the figure, the element values are given 
as percents. The elements covered by the Timing Advance zone are given, 
the value of 100, while the other elements are given the value of zero. 

FIG. 5 illustrates the estimation process performed by the accu- 
5 racy server. As mentioned above, upon receiving a parameter set the accu- 
racy server calculates a matrix corresponding to this parameter set or re- 
trieves said matrix from the database 39a using the parameter set as a 
search key (step 51). In the latter case the matrices are preferably formed 
using a suitable calibration mechanism, such as the GPS system. In other 

10 words, when the matrices are formed before the commissioning of the sys- 
tem, the distribution of the location of the mobile is determined by the calibra- 
tion mechanism and the matrix is formed on the basis of the results given by 
the calibration mechanism. A combinatory mechanism can also be used, for 
example, so that the matrices calculated in advance are further processed in 

15 the accuracy server in order to obtain the final matrix describing the probabil- 
ity distribution. In a preferred embodiment of the invention, this is done by 
using map information to adjust the element values of the matrix. 

In this embodiment the accuracy server retrieves map information 
in order to weight the element values by this data. The matrices to be 

20 weighted can originate from a database or can be formed in real-time by the 
accuracy server. Map information refers here to any information which corre- 
lates with the location probability or with the movement probability of the 
mobile. The map information used by the accuracy server is preferably in the 
form of map matrices similar to the matrices formed on the basis of the pa- 

25 rameter sets, except that in the map matrices an individual element indicates 
the type of the surface covered by the element. The classification into surface 
types can be simple, such as road, forest, river, etc., each type being as- 
signed a weight value. Some types are stronger than others; for instance, if a 
road goes through a forest, the subscriber is more likely traveling along the 

30 road rather than through the forest, so the type of the element is set to "road". 
The map matrix can be formed in real-time from the information retrieved 
from the map database, or the information in the database can be in the form 
of matrices. 

The element values of the matrix formed on the basis of the pa- 
35 rameter set received are then weighted according to the type of land in ques- 
tion (step 53), whereby a weighted matrix is obtained. A previous matrix 
obtained in connection with a preceding parameter set relating to the same 
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mobile is also retrieved from the history database 39b, and this matrix is 
supplied to an update process. In this context, the previous matrix is also 
called the history matrix. In the update process (step 54) the effect of the 
estimated movement of the mobile on the location distribution is taken into 
5 account by spreading the probability distribution of the previous matrix ac- 
cording to the movement data. This operation is based on the information 
available about the movement of the mobile, such as speed and/or direction 
information. Map correlation coupled with movement information can be used 
to yield more accurate results. If no information on the movement is avail- 

10 able, the probability distribution is spread evenly in all directions. As a result 
of the update process, a spread history matrix is obtained. 

The current matrix and the spread history matrix are supplied to a 
combination process where a combined matrix of the same size is deter- 
mined according to predetermined rules (step 55). At least one location esti- 

15 mate is then calculated on the basis of the combined matrix (step 56). Fur- 
thermore, the combined matrix is supplied to the history database, where it 
replaces the previous matrix, i.e. the next time the previous matrix is re- 
trieved, the combined matrix calculated in connection with the then preceding 
parameter set is obtained as the previous matrix to be updated. 

20 FIG. 6 illustrates the updating of the previous matrix, in which the 

values of the previous matrix are updated according to the assumed move- 
ment of the mobile. The update process includes one or more calculation 
rounds. During each round the probability distribution carried by the history 
matrix is spread once. The process employs three parameters which define 

25 how many rounds are needed each time; n, T, and T1 , where n indicates the 
number of times the spreading has been performed, T indicates the time 
elapsed from the previous update process (i.e. the time elapsed since the 
location of the mobile was estimated previously), and T1 is a fixed time inter- 
val which is deducted from T every time the matrix is spread. 

30 The process first forms n1xn1 basic coefficients Wg used in the 

spreading process (step 60). In this example n1=3 and all the coefficients get 
the value of 1 as follows: 
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Parameter n is then given the value of zero (step 61), since no 
spreading has yet taken place, and the value of T1 is deducted from the 
value of T (step 62). 

The basic coefficients are then adjusted according to the move- 
5 ment information relating to the mobile, such as the direction and speed 
information (step 63). Map information can also be utilized in determining the 
adjusted coefficients. It is assumed in this example that the mobile is moving 
directly eastwards, whereby the adjusted coefficients w'y are as follows: 
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10 In a spreading process the adjusted coefficients are then used to 

adjust the value of each element e K j of the history matrix as follows (step 64): 



where E kJ is the element value of the spread matrix and M is the 
sum of adjusted coefficients deviating from zero. Thus the spread matrix is 
15 obtained so that each element of the previous matrix is replaced by a 
weighted average of itself and those of its neighbors which do not get a zero 
weight. In practice the number of elements affecting the adjusted value of 
each element can be greater, for example 16 (4x4 weights) or 25 (5x5 
weights). 

20 After the first spreading round, the accuracy server tests whether 

T is still greater than zero and n smaller than a predefined limit n max (step 65). 
If this is the case, the above-described spreading steps 63 and 64 are re- 
peated. At the beginning of each new round, the value of n is incremented by 
one (step 66), and the value of T1 is deducted from the current value of T 

25 (step 62). A new matrix is calculated as long as T remains greater than zero, 
and n smaller than a predefined limit n max . 

As to step 63, it is further to be noted that the adjustment of the 
basic coefficients by means of the movement information is performed only 
once for each update process (i.e once for each step 54). The adjustment 

30 according to the map information has to be made separately for each matrix 
element during each spreading round, since the map information used to 
adjust the coefficients is element-specific, i.e. depends on the area defined 
by the element. 
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FIG. 7 illustrates how the matrix of FIG. 4 has changed after one 
spreading round when the above basic coefficients are used. Thus in this, 
case no direction-based weighting is employed, but it is assumed that all 
directions are equally probable. As can be seen from the figure, the probabil- 
5 ity distribution has spread; the number of elements having a non-zero value 
has increased from 17 to 45. 

FIG. 8 illustrates how the matrix of FIG. 4 has changed after one 
spreading round when the above-described adjusted coefficients are used. 
Thus, in this case it is assumed that the mobile is traveling eastwards. Now 

10 the number of elements having a non-zero value has increased to 37. Utiliz- 
ing the map information in the spreading process means that the probabilities 
spread faster along roads than along forests, for example. 

The spread history matrix and the current matrix can be combined 
in various ways. However, there are always two probability values, A and B, 

15 corresponding to a certain area, one from the element corresponding to the 
respective area in the current matrix and another from the element corre- 
sponding to the respective area in the spread matrix. In the combination 
process (step 55), a new element value is determined for the respective 
element of the combined matrix based on these two values. There are at 

20 least three possible ways to determine the new element value for each ele- 
ment of the combined matrix: 

1. calculate the sum of the two values, i.e. A+B, 

2. calculate a weighted sum of the two values, i.e. rxA+(1-r)xB, 

3. select the value which is the highest one, i.e. max{A, B}. 

25 As mentioned above, the combined matrix formed in one of the 

above ways is stored as a new history matrix. 

Based on the values of the combined matrix, the final estimate(s) 
can also be determined in various ways. At least the following ways are pos- 
sible: 

30 1. The coordinates of the element containing the highest value are 

selected as the location estimate. 

2. The coordinates corresponding to the center of gravity of the 
element values are selected as the location estimate. The center of gravity 
here refers to an average (calculated in a desired way) of the coordinate 

35 values associated with the element values. 

3. A certain area, which is obtained based on the previous esti- 
mate and the speed and/or direction information, is examined. The element 
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containing the highest value within that area is then selected as the location 
estimate. Another alternative is to select as the location estimate the coordi- 
nates corresponding to the center of gravity of the element values of the 
area. 

5 4. Local averages can be calculated for element groups of the de- 

sired size. For example, averages can be computed over local areas of 
100x100 m 2 , and the center point of the local area having the highest aver- 
age can be selected as the local estimate. The form of the local area can 
also be other than a square, for example a circle. 
10 5. Areas of predetermined form, such as circles, are drawn around 

estimate candidates, each area covering a certain proportion of the total sum 
of the elements. The candidate is then selected for which the said area is the 
smallest. 

The best alternative can be selected by testing the system in the 

15 desired environment and finding the alternative yielding the best results. 

The formation of the matrix corresponding to the current parame- 
ter set (step 51) is discussed next. As mentioned above, the matrix can be 
formed directly in real-time on the basis of the parameter set, or it can be 
retrieved from among the matrices produced in advance by using the pa- 

20 rameter set as a search key. A third alternative is to use a combination of 
these two methods. 

When the matrix is formed in real-time from the parameter set, the 
geometrical counterpart of the parameter set, i.e. the area defined by the 
parameter set, can be superimposed on an empty matrix covering the area in 

25 question. FIG. 9 illustrates this by showing a Timing Advance zone 90 super- 
imposed on a 7x7 matrix. Each of the elements is then assigned a value 
which is directly proportional to the proportion which the said counterpart 
covers of the area of the element. In the example of FIG. 9 the values range 
from one (no coverage) to ten (full coverage). These values can then be 

30 adjusted employing the map information. Alternatively, all the elements over- 
lapping the geometric counterpart can be assigned the value of one, while 
the rest of the elements, which do not overlap with the element, are assigned 
the value of zero. 

When the matrix is retrieved from a database containing matrices 

35 formed in advance, it is preferable to measure the actual probability distribu- 
tion by means of a suitable mechanism, such as the GPS. Thus, in this case 
the matrix elements indicate the probability distribution measured in advance 
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in real environment. The parameter sets which are substantially equal to 
each other, i.e. which with a desired accuracy have the same parameter 
values, are associated with the same matrix. 

The third alternative which combines the above mechanisms could 
be such that an in-advance calculated matrix is always used when it is avail- 
able, while the matrix is computed in real-time only when no calculated matrix 
is available in advance. This is because a matrix determined in-advance 
probably represents a more reliable perception of the probability distribution. 

Although the invention was described above with reference to the 
examples shown in the appended drawings, it is obvious that the invention is 
not limited to these, but may be modified by those skilled in the art without 
departing from the scope and spirit of the invention. In a simplified form of the 
invention, no map information or direction information is used, i.e. the prob- 
ability distribution is spread evenly in all directions. More than one history 
matrix can be employed for processing the current matrix. The network ele- 
ments used can also vary. For example, the parameter sets can be received 
from any network entity having access to them, and the steps of the method 
can be distributed among different network elements. In the future a mobile 
may perform some or all of the steps of the method. Depending on the con- 
tent of the parameter set, a mobile might simply need external information 
required to form the matrices in order to be able to perform the steps of the 
method. . - 
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Claims 

1. A method for locating mobile terminals in a mobile network, the 
method comprising the steps of: 

receiving location-dependent parameter sets, each parameter set 
5 to include at least one parameter indicative of the location of an individual 
mobile terminal, and 

determining a location estimate for a parameter set received, the 
location estimate indicating the location of the respective mobile terminal, 
characterized by the steps of: 
10 for an individual parameter set, forming a matrix corresponding to 

said set, the matrix comprising a plurality of elements, whereby each element 
is associated with a certain geographical area and contains a value indicating 
a probability of the mobile being located within said area, 

storing at least one matrix formed for a mobile, and 
15 in response to a parameter set received for the mobile, 

- updating the values of at least one matrix stored for the 

mobile and 

- determining the location estimate on the basis of the ele- 
ment values of the matrix corresponding to the parameter set received and 

20 on the basis of the element values of the said at least one matrix having the 
updated values. 

2. The method as defined in claim 1, characterized in that 
the determining step includes 

- combining the element values of the matrix corresponding to the 
25 parameter set received and the element values of said at least one matrix 

having the updated values according to predetermined rules, whereby a 
combined matrix is obtained and 

- defining the location estimate on the basis of the combined ma- 
trix. 

30 3. The method as defined in claim 2, characterized in that 

the storing step includes storing the combined matrix and 
the updating step includes updating the values of the combined 

matrix obtained in connection with a preceding parameter set received for the 

mobile. 

35 4. The method as defined in claim 3, characterized in that 

the updating step includes updating the values repetitively in successive 
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calculation rounds, whereby the updated values obtained in a calculation 
round are updated in the next calculation round. 

5. The method as defined in claim 4, characterized in that 
the number of calculation rounds is directly proportional to the time elapsed 

5 since the updating step performed in connection with the preceding parame- 
ter set received for the mobile. 

6. The method as defined in claim ^characterized in that 
the forming step includes weighting said values on the basis of map informa- 
tion describing surface types of said geographical areas. 

10 7. The method as defined in claim 1, characterized in that 

the updating step includes weighting the values of said at least one matrix on 
the basis of information indicating the movement of the mobile. 

8. The method as defined in claim 1, characterized in that 
the updating step includes weighting the values of said at least one matrix on 

15 the basis of map information describing surface types of said geographical 
areas. 

9. The method as defined in claim 1, characterized in that 
the forming step includes calculating the matrix in response to a parameter 
set received. 

20 10. The method as defined in claim 1, characterized in that 

the forming step includes 

- calculating a plurality of matrices in advance, 

- associating each matrix with a parameter set, and 

in response to a parameter set received, retrieving a matrix corre- 
25 sponding to said parameter set. 

1 1. A system for locating mobile terminals in a mobile network, the 
system comprising: 

first means for receiving parameter sets, each parameter set in- 
cluding at least one parameter indicative of the location of an individual mo- 
30 bile terminal and 

second means for finding a location estimate for a parameter set 
received, the location estimate indicating the location of the respective mobile 
terminal, 

characterized in that the system comprises 
35 third means for forming a matrix corresponding to a parameter set, 

the matrix comprising a plurality of elements, whereby each element is asso- 
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ciated with a certain geographical area and contains a value indicating a 
probability of the mobile being located within said area, and 

fourth means for storing at least one matrix formed for a mobile, 
wherein the second means, responsive to a parameter set re- 
5 ceived for the mobile, are adapted (a) to update the values of at least one 
matrix stored for the mobile and (b) to determine the location estimate on the 
basis of the element values of the matrix corresponding to the parameter set 
received and on the basis of the element values of the matrix with the up- 
dated values. 

10 12. A mobile terminal for a mobile network, the mobile terminal 

comprising: 

first means for receiving parameter sets, each parameter set in- 
cluding at least one parameter whose value is dependent on the location of 
the mobile terminal within the network and 
15 second means for finding a location estimate for a parameter set 

received, the location estimate indicating the location of the mobile terminal, 
characterized in that the mobile comprises 
third means for forming a matrix corresponding to a parameter set, 
the matrix comprising a plurality of elements, whereby each element is asso- 
20 ciated with a certain geographical area and contains a value which indicates 
a probability of the mobile locating within said area, and 

fourth means for storing at least one matrix in the mobile, 
wherein the second means, responsive to a parameter set re- 
ceived, are adapted (a) to update the values of at least one matrix stored and 
25 (b) to determine the location estimate on the basis of the element values of 
the matrix corresponding to the parameter set received and on the basis of 
the element values of the matrix with the updated values. 

13. A computer program product stored on a computer readable 
storage media, the product being adapted to perform the steps of claim 1 
30 when run on a computer. 
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Abstract 

The invention concerns locating mobile terminals in a mobile network. A 
location estimate is determined based on a parameter set received from the 
mobile network. In order to improve the accuracy of the system, a matrix is 
5 formed for an individual parameter set, the matrix comprising a plurality of 
elements, whereby each element is associated with a certain geographical 
area and contains a value which indicates the probability of the mobile being 
located within said area. At least one matrix formed for a mobile is stored, 
and in response to a parameter set received for the mobile, (a) the values of 
10 at least one matrix stored are updated, and (b) the location estimate is de- 
termined on the basis of the element values of the matrix corresponding to 
the parameter set received and on the basis of the element values of the said 
at least one matrix having the updated values. 
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